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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


TECHNICAL  NOTE  2605 

BEHAVIOR  OF  VORTEX  SYSTEM  BEHIND  CRUCIFORM  WINGS  - 
MOTIONS  OF  FULLY  ROLLED -UP  VORTICES 


By  Alvin  H.  Sacks 


SUMMARY 


The  motions  of  four  fully  rolled-up  vortices  representing  the  vortex 
system  trailing  behind  cruciform  wings  are  studied  by  theoretical  and 
visual-flow  methods.  The  analysis  applies  throughout  the  Mach  number 
range  * 


Equations  are  developed  for  the  three-dimensional  paths  traced  by 
the  vortices  behind  a  cruciform  wing  banked  ^5° ,  and  calculations  are 
made  of  the  distance  behind  the  wing  at  which  the  upper  two  vortices 
pass  through  the  lower  two.  It  is  found  that  this  "leapfrog"  distance 
depends  upon  the  lift  coefficient,  aspect  ratio,  and  span  loading  of  the 
cruciform  wing,  and  that  for  low-aspect -ratio  cruciform  wings  leapfrog¬ 
ging  may  occur  within  two  chord  lengths  of  the  trailing  edges. 

The  various  types  of  vortex  motion  to  be  expected  throughout  the 
angle -of -attack  range  are  considered  in  some  detail,  and  the  interaction 
of  the  two  vortex  sheets  shed  from  the  cruciform  wing  is  taken  into 
account.  Results  of  some  water-tank  studies  are  also  presented  and 
compared  with  the  theory. 


INTRODUCTION 


The  downwash  behind  plane  wings  has  been  studied  theoretically  by  a 
number  of  authors  and  considerable  attention  has  been  given  to  the 
rolling  up  of  the  trailing  vortex  sheet.  The  analysis  of  Kaden 
(reference  1)  predicted  the  distance  behind  the  wing  at  which  the  sheet 
may  be  considered  to  be  fully  rolled  up  into  two  trailing  vortices,  and 
this  work  was  later  used  (in  reference  2)  to  demonstrate  the  usefulness 
of  the  single  horseshoe-vortex  approximation  for  the  calculation  of  the 
downwash  behind  wings  of  low  aspect  ratio.  Since  at  the  present  time 
cruciform  configurations  are  largely  confined  to  wings  of  low  aspect 
ratio,  the  rolling  up  of  the  trailing  vortex  sheets  is  again  of  major 
concern,  and  the  behavior  of  the  fully  rolled-up  vortices  is  again  of 
considerable  interest. 
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While  in  the  case  of  the  plane  wing  the  vortex  sheet  became  rolled 
up  into  a  vortex  pair  which  simply  moved  downward  at  a  uniform  speed  at 
great  distances  behind  the  wing,  the  analogous  problem  for  cruciform 
wings  is  necessarily  more  complicated.  Instead  of  the  two  rolled-up 
vortices,  there  are  now  presumably  four  (one  from  each  wing  panel)  and 
their  induced  effects  upon  one  another  are  such  as  to  produce  quite 
intricate  paths  of  motion.  The  downwash  field,  of  course,  may  therefore 
become  extremely  involved.  This  report  is  concerned  with  a  study  of  the 
motions  of  the  four  rolled-up  vortices  and  their  effect  on  the  downwash 
behind  cruciform  wings. 
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SYMBOLS 


aspect  ratio 
span  of  one  wing  (2s0) 
root  chord 

lift  coefficient  of  cruciform  wing 

lift  coefficient  of  plane  wing  (  — 

\  qS 

section  lift  coefficient  for  plane  wing 
distance  behind  wing  trailing  edge 

distance  behind  wing  trailing  edge  where  four  rolled-up  vortices 
are  collinear  (leapfrog  distance) 

constant  of  the  motion  related  to  distance  between  centers  of 
gravity  of  two  vortex  systems  (yx  +  y2) 

parameter  which  depends  on  initial  vortex  positions 


lift  of  cruciform  wing 
lift  of  plane  wing 
free -stream  Mach  number 


free -stream  dynamic  pressure  plT* ^ 

parameter  used  for  periodic  motion  — — - 
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Q2  parameter  used  for  aperiodic  motion  (-Q2) 

t 

R2 

4(g  +  4) 

S  area  of  one  wing 

'  s  local  semi span  of  one  wing 

sQ  maximum  semi  span  of  one  wing 

sT  semispan  of  rolled-up  vortices  from  one  wing 

t  time 

U  free-stream  velocity 

v  velocity  of  vortex  center  in  y  direction 

w  velocity  of  vortex  center  in  z  direction 

x,y, z  right-hand  orthogonal  coordinates  with  x  axis  in  stream 
direction 

^  y0,  zQ  initial  values  of  y1  and  z± 

zi>2  vertical  distance  between  vortices  1  and  2 
(See  sketch  on  page  9.) 

a  angle  of  attack  of  cruciform  wing  axis 

aT  angle  of  attack  of  one  wing 

|3  angle  of  sideslip  of  one  wing 

7  ratio  of  vortex  strengths 

r  circulation  (positive  counterclockwise) 

p  mass  density 

Subscripts  1,  2,  3,  and  4  refer  to  vortex  number.  (See  sketch  on  page  5. 
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GENERAL  ANALYSIS 


The  phenomenon  of  the  rolling  up  of  the  vortex  sheet  behind  plane 
wings  was  discussed  as  early  as  1907  by  Lanchester  (reference  3)  and  has 
since  been  the  subject  of  a  large  number  of  investigations.  As  a  result 
of  such  investigations,  it  has  been  found  (reference  2)  that,  whereas 
the  vortex  sheet  behind  wings  of  high  aspect  ratio  may  be  considered  to 
remain  flat  at  the  tail  location,  a  good  approximation  to  the  downwash 
several  chord  lengths  behind  wings  of  low  aspect  ratio  is  often  obtained 
by  considering  the  sheet  to  be  fully  rolled  up  into  two  trailing  vortices. 
In  the  case  of  the  cruciform  wing,  as  in  the  case  of  the  plane  wing,  the 
rate  of  rolling  up  of  the  vortex  sheets  will  depend  upon  the  aspect  ratio 
of  each  wing.  It  should  therefore  be  possible  to  obtain  a  reasonably 
good  approximation  to  the  downwash  field  behind  low-aspect-ratio  cruci¬ 
form  wings  by  considering  the  vortex  sheets  to  be  fully  rolled  up;  that 
is,  by  replacing  the  two  vortex  sheets  by  four  discrete  vortices.  It  is 
with  this  simplified  model  of  the  physical  problem  that  the  analysis  of 
this  report  is  concerned  at  the  outset.  Subsequently,  it  will  be  seen 
that  the  theory  can  be  modified  to  provide  a  more  accurate  representation 
of  the  actual  flow  field. 

With  low-aspect-ratio  wings  in  source-free  flow  the  linearized  dif¬ 
ferential  equation  for  the  perturbation  velocity  potential  9 


(l-M2)9  +9  +9  =  0 

'  /  1  w  1  inr  rr  rr 


XX 


yy 


zz 


(1) 


can  be  satisfactorily  approximated  at  all  Mach  numbers  by  the  two- 
dimensional  equation 

<Pyy  +  <PZZ  =  0  (2) 

since  (l-M2)  cp^  is  much  smaller  than  Tyy  and  <PZZ  if  o/l-M2  A  is 
small.  (See  references  k  and  5.)  Therefore,  the  problem  of  the  motions 
of  the  four  rolled-up  vortices  can  first  be  treated  by  latera.1  strip 
theory  as  a  two-dimensional  problem  in  planes  perpendicular  to  the 
flight  direction.  For  the  cruciform  wing  at  ^5°  angle  of  bank,  where  a 
vertical  plane  of  symmetry  is  present,  the  problem  is  thus  reduced  to 
that  of  the  motion  of  two  pairs  of  vortices  with  a  common  axis.  This 
problem  was  first  attacked  by  Love  (reference  6)  and  later  by  Hicks 
(reference  7)>  both  in  connection  with  the  analogous  three-dimensional 
problem  of  the  motion  of  vortex  rings.  Unfortunately,  neither  of  these 
authors  worked  completely  the  problem  of  interest  here,  the  former  being 
concerned  only  with  the  relative  paths  of  the  vortex  pairs  about  each 
other,  whereas  the  latter  was  concerned  chiefly  with  the  actual  problem 
of  the  vortex  rings. 


Since  this  paper  is  concerned  with  the  three-dimensional  paths  of 
the  four  rolled-up  vortices  behind  the  wing,  the  paths  in  transverse 
planes  will  first  be  determined,  and  will  then  be  related  to  the  elapsed 
time  or  the  distance  the  wing  has  traveled  in  the  flight  direction. 

The  analysis  will  be  confined  to  the  case  of  ^5°  bank. 
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Determination  of  Vortex  Paths 


Relative  motions  in  transverse  plane.-  The  motion  of  any  system  of 
two-dimensional  vortex  filaments  is  determined  by  the  induced  velocity  of 
each  vortex  due  to  all  the  other  vortices  of  the  system.  It  is  shown  in 
reference  8  that  the  induced  velocity  of  any  particular  filament  can  be 
obtained  from  the  function 


W  =  ^  Ermrn  log  ( (3) 


where  m<n  and  =  Ym  +  izm.  That  is,  if  one  writes 

-W  =  ®  +  it  (4) 


then  ^  is  a  constant  of  the  motion  and  is  analogous  to  the  stream 
function  in  giving  the  components  of  the  velocity  of  the  particular 
vortex.  Thus,  the  path  of  each  vortex  can  "be  determined  from  equa¬ 
tions  (3)  and  (4),  provided  that  the  motion  of  one  vortex  determines  the 
motions  of  all  the  other  vortices  of  the  system.  (Such  is  the  case,  for 
instance,  where  all  the  vortices  but  one  are  images  of  the  one.)  In  the 
present  problem,  since  only  two  of  the  four  vortices  are  images,  equa¬ 
tions  (3)  and  (4)  are  not  sufficient  to  describe  completely  the  motion  of 
the  system.  However,  these  equations  will  be  useful  in  determining  the 
paths  of  two  of  the  vortices  relative  to  the  other  two.  It  will  be  seen 
later  that  this  information  is  of  considerable  interest. 


For  the  cruciform  wing  banked  45°,  the  simplified  vortex  system  as 
viewed  in  the  y,z  plane  consists  of  two  pairs  of  two-dimensional 
vortices  having  a  common  plane  of  symmetry  as  shown  in  the  sketch. 

It  is  apparent  from  the  symmetry  of 
the  configuration  that 


ri  -  -r4; 


r2  =  -r3 


Further,  if  each  wing  of  the  cruciform 


3 


can  be  considered  as  a  plane  wing  in 

sideslip  (as  was  demonstrated  in  I 
reference  9  for  the  calculation  of  | 
load  distribution)  it  follows  that  [ 


since  the  vortex  system  far  behind  a 
plane  lifting  wing  must  consist  of  two 
equal  (and  opposite)  vortices, 
although  their  positions  will  depend 
upon  the  angle  of  sideslip  while 
their  strength  will  depend  upon  the 


- — — — - ►  y 

plane  of 
symmetry 
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lift.  On  this  basis,  it  can  be  concluded  that 

rx  =  r2  =  -r3  =  -r4  =  r 

Actually,  since  the  vortex  sheets  trailing  behind  the  cruciform  wing 
must  influence  one  another  during  the  rolling-up  process  itself,  the 
assumption  that  the  wings  act  independently  is  not  strictly  applicable 
to  considerations  of  the  wake  as  in  the  present  problem.  However,  this 
complication  will  be  deferred  to  a  later  section  and  the  analysis  will 
first  be  developed  on  the  basis  of  four  vortices  of  equal  strength. 
Hence,  from  this  point  on,  the  subscripts  will  be  omitted  in  reference 
to  the  vortex  strengths  which  will  be  referred  to  simply  as  plus  or 
minus  T . 

When  the  directions  of  rotation  of  the  four  vortices  are  taken  into 
consideration,  equation  (3)  can  be  written  for  the  present  problem  as: 

W  =  7I3  [log  (^1-^2)  -  log  (£2-^3)  +  log  (^3-U)  - 


log  (£i-£4)  -  log  ( £i-£3)  -  log  (£2-C4)  (5) 


Noting  further  that  £4  and  £3  are  simply  the_negative  complex  conju¬ 
gates  of  and  £a,  that  is  £4  =  -£x,  £3  =  equation  (5)  can  be 

expressed  as 

iP2  f  _  _ 

■  w  =  log  (Ci-(;2)  -  log  (C2  +  ls)  +  log 

log  (Ci  +  lx)  -  log  (Cx  +  Is)  -  log  (£a  +  lx) 

-i£2iog  ikMkk  1 

^  L(C2+^2)(^+Ci)(Cx+C2)(52+Ci)  J  (6) 

Due  to  the  symmetry  (see  sketch  on  page  5)^  the  vortex  system  may 
also  be  considered  as  two  groups  of  vortices  with  equal  and  opposite  total 
circulation  (one  group  on  either  side  of  the  plane  of  symmetry).  For 
such  a  pair  of  groups,  the  centers  of  gravity1  of  the  two  groups  must 
remain  a  fixed  distance  apart.  (See  reference  10.)  In  other  words, 
the  center  of  gravity  of  each  group  must  move  parallel  to  the  plane  of 
symmetry.  Since  here  the  strengths  are  equal  in  magnitude,  this  may  be 
stated  as 

zyin  r(yx+y2)  yi+y2  f  _  _+QT1+ 

yc.g.  =  -sir"  — 2T  =  2  s  5 "  constant  w) 


■^The  center  of  gravity  of  a  group  of  vortices  is  defined  as  the  center 
of  gravity  of  a  similar  field  of  point  masses,  the  mass  of  each  being 
proportional  to  the  strength  of  the  corresponding  vortex.  Negative 
masses  correspond  to  negative  circulations.  _ 
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With  this  information,  equation  (6)  can  he  expressed  in  terms  of  the 
coordinates  y  and  z  and  simplified  to 


Now,  since  W 


w  =  •*Jg  log  I (yi-yg)  +  (zi~za)g  \ 

1  ViJa .  [  f*+ ( z1  -z£  )2]  J 


is  purely  imaginary,  equation  (4)  is  simply 

(yi-ya)g  +  (zx-z g)g  1 
^yiyet^+Czi-z 2)23  J 


(8) 

(9) 


At  this  point,  a  moving  coordinate  system  is  introduced  such  that  one 
pair  of  vortices  always  lies  on  the  y  axis.  This  is  done  hy  a  change 
of  variable 

Zl,2  "  zi  "  za  (10) 

Now,  noting  that  y2  =  f  -  yx,  equation  (9)  becomes 


(2y2-f)a  +  Zl.aa  ' 

%1(f-y1)(f2+z^a2)  _ 


(11) 


Setting  'If  of  equation  (ll)  equal  to  a  constant  then  yields  an  equation 
for  the  paths  of  vortices  1  and  2  relative  to  each  other#  Thus, 


(2yi-f)g  +  zi,gg 
y^f-yj^+z^/) 


=  constant 


or 


determines  the  relative  paths  of  vortices  1  and  2  (and  also  3  and  %)  once 
the  value  of  the  constant  G  is  determined . 


In  order  to  evaluate  G  of  equation  (12),  one  must  determine  the 
initial  positions  of  the  four  vortices.  It  must  he  realized,  of  course, 
that  initially  (i.e.,  at  the  trailing  edges  of  the  cruciform  wings)  there 
are  not  four '‘discrete  vortices  hut  a  cruciform  system  of  vortex  sheets. 
The  immediate  problem,  then,  is  to  replace  these  sheets  by  four  vortices 
properly  located  in  the  £(or  y,z)  plane.  For  the  present,  it  will 
suffice  to  observe  that  the  initial  locations  of  the  four  vortices  must 
lie  on  the  straight  lines  formed  by  the  trailing  edges  of  the  wings. 
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Thus,  for  45°  hank, 


zi ,2  =  yi  +  y2  =  f  When  yx  =  yQ  (13) 

where  the  subscript  o  refers  to  the  initial  value  of  y1 .  The  actual 
determination  of  y0  will  he  deferred  to  a  later  section  of  this  report 
in  which  the  physical  problem  will  be  considered  in  more  detail. 

With  the  boundary  condition  expressed  in  equation  (13),  equa¬ 
tion  (12)  can  be  written  in  terms  of  the  initial  value  of  yx : 


It  will  be  shown  later  that  y0,  and  hence  G,  is  a  function  of  angle 
of  attack  for  the  cruciform-wing  configuration. 

While  equation  ( ih-)  does  not  completely  describe  the  motion  of  the 
four  vortices,  it  yields  some  interesting  information  regarding  the  types 
of  motion  to  be  expected.  For  instance,  if  the  relative  path  described 
by  equation  (l4)  is  a  closed  curve  it  can  be  concluded  that  the  motion 
is  periodic.  The  relative  paths  of  the  vortices  are  plotted  in 
figure  1  for  various  values  of  G,  and  it  is  found  that  the  motion  is 
periodic  for  G  less  than  4  and  aperiodic  for  G  greater  than  4. 

This  is  in  agreement  with  the  findings  of  Love  (reference  6)  who 
showed  that  the  condition  for  periodic  motion  is  that 

r  <  3  +  2  */2~  (15) 

where  r  =  ~  at  z,  „  =  0.  That  is,  r  is  the  ratio  of  the  lateral 
Yl  1,2 

displacements  from  the  plane  of  symmetry  when  the  four  vortices  are 
collinear  (when  one  pair  is  passing  through  the  other) .  This  type  of 
vortex  motion  was  first  discussed  "by  Helmholtz  in  connection  with 
vortex  rings  having  the  same  axis  and  circulations  in  the  same 
direction  (reference  11): 

"The  foremost  widens  and  travels  more  slowly,  the  pursuer 
shrinks  and  travels  faster,  till  finally,  if  th^ir  veloc¬ 
ities  are  not  too  different,  it  overtakes  the  first  and 
penetrates  it.  Then  the  same  game  goes  on  in  the  opposite 
order,  so  that  the  rings  pass  through  each  other  alternately 
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Complete  motions  in  transverse  plane .  -  Since  the  relative  paths  of 
the  vortices  on  either  side  of  the  plane  of  symmetry  are  given  by  equa¬ 
tion  (l4),  establishing  the  actual  path  of  any  one  of  the  four  vortices 
will  completely  determine  the  motion  of  the  system  in  the  y, z  plane. 
Therefore,  the  actual  path  of  vortex  1  will  now  be  calculated. 


The  determination  of  any  of  the  individual  vortex  paths  requires 
the  solution  of  the  following  equation: 


Zn 


dyn/dt 


(16) 


where  vn  and  wn  are  the  velocity 
components  of  the  nth  vortex  in  the 
y  and  z  directions,  respectively. 
If  equation  (l6)  is  evaluated  for 
one  vortex,  say  vortex  1,  then  the 
path  of  vortex  2  can  be  obtained 
from  equations  (7)>  (10),  and  (lb), 
and  the  paths  of  vortices  3  and  k 
are  found  by  symmetry.  It  will  be 
seen  that  this  is  considerably 
easier  than  solving  four  equations 
of  the  type  given  by  equation  (l6) . 
By  adding  the  contributions  of  all 
the  other  vortices  to  the  velocity 
components  of  vortex  1,  one  finds 
( see  sketch) : 


and 


Vj_  =  -  -  COS  0^  +  -  COS  @3 

2Tcr2  2*r3 

r  r  r 

w«l  as  -  -  -  — —  sin  @2  -  -  sin  03  (17) 

2irr4  2jtr2  2jtr3 


Further,  from  the  geometry  of  the  system,  it  is  seen  that 


COS  ©2 

Zl,2 

COS  03 

Zl,2 

r2 

'  zi,22  +  (f  -  2yi)2' 

r3 

Zl,22  +  f2 

sin  02 

f  -  2yi 

sin  ©3 

f 

r2 

zl,  22  +  (f  -  2Yl)2’ 

r3 

Zl,22  +  ^ 

and 
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The  velocity  components  va  and  vx  can  now  be  expressed  in  the  following 
form: 

2rzi,2  r  yi(f-yi)  1 


vi  =  - 


it  \[Zl>22  +  (f-2yi)2][Zl>22  +  f2] 
T  f  (f2  +  2i,gg)g-  8fyx2  (f  -  yj 

1  ly^2  +  zi,22)tzi,22  +  (f  ~  2yx )21 


If  equations  (l4)  and  (18)  are  used  to  express  vx  and  Wj  as  functions 
of  yx ,  then  equation  (l6)  can  be  integrated  directly  to  give  the  path 
of  vortex  1.  After  algebraic  manipulation,  equation  (16)  thus  becomes 


j1  f4+2f3y1(G-l)-Gfgyi2(G+2)+2Ggfy13-Ggy14 _ 

y0  (fg-Gfy1+Gy12) */[ -(G+4)y12+f (G+4)y1-f2] [ Gyx2 -Gf yx  +f2 ] 


Zi-Z.  = 


where  zQ  is  the  initial  value  of  z1. 

If  now  the  numerator  of  the  integrand  is  divided  by  the  rational 
factor  in  the  denominator,  and  the  substitution 

*  -  ♦  i 

is  made,  then  the  quantity  (zx  -  z0)  can  be  expressed  as  the  sum  of  four 
integrals 


Z  i  "’Z  — 


yG(G+4) 


y  -£ 
•'i  2 


-  Q2G 


4  V(R2-ti2)(qW) 


.I  */  (R2-t]2)  (Q2+t)2) 


f 

.4  Ji-5 


f 

2f 3  o1'2 


GJy  J  (Q2+q2)  <y(R2-Tl2)(Q2+n  )  G  y 

^02  jO  2 


2  f%  2  f2(4-G) 

where  R  =  -  and  Q  =  - 

4(g+4)  4g 
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Equation  (20)  is  now  in  a  convenient  form  if  G  is  less  than  4.  This 
corresponds  to  the  range  for  periodic  motion.  On  the  other  hand,  if  G 
is  greater  than  4  (aperiodic),  a  new  constant 

4*  -  -Q2  = 

is  introduced  so  that  equation  (20)  becomes 


third  possibility  is  that  G  is  equal  to  4,  in  which  case 
Q  =  Q2  =  0.  Equations  (20)  and  (21)  then  both  reduce  to 


The  three  cases  corresponding  to  equations  (20),  (21),  and  (22)  will 
be  treated  separately  since  the  integrations  will  be  somewhat  different 
for  each. 

For  cases  where  the  vortex  motion  is  periodic,  that  is,  when  G  is 
less  than  4,  the  determination  of  the  vortex  paths  requires  the  evalua¬ 
tion  of  the  integrals  of  equation  (20) .  It  is  seen  that  the  fourth 
o-f  this  equation  is  elementary  while  ^the  first  three  are 
•  Thus,  while  the  fourth  integral  can  be  evaluated  simply  by 

noting  that  q  dt|  =  1  d(qe),  the  other  integrals  require  the  use  of  the 

elliptic  transformations2 


'These  and  subsequent  transformations  were  obtained  from  a  comprehensive 
table  of  elliptic  integrals  prepared  by  Paul  Byrd,  Ames  Laboratory, 
MCA. 
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sn2  ui  =  sin2  *Pi  = 


R2- v  yi- 


sn2  uQ  =  sin2  90 


RC“  W  I 


P  P  ✓*  '  ' 

r2  +  or  16 

With  these  transformations  the  integrals  can  be  evaluated  and  the -results 
expressed  in  terms  of  elliptic  integrals  of  the  first  and  second  kinds. 
The  final  solution  of  equation  (20),  after  collection  of  terms  and  eval¬ 
uation  of  the  constants,  is  given  in  nondimensional  form  by  the  expres- 


E(k,  T0)-E(k,  Ti) 


sin  9i  cos  cpi  sin  To  cos  To 


8(4-G)  Wl-k  sin2  Ti  vl-k2  sin2  Tc 


g  a/g(g+4) 


sin  Ti 


1-k2  sin2  Ti 


sin  Tr 


l-k2sin2  T, 


•where  k  =  and  it  is  recalled  that  zQ  is  equal  to  yQ  due  to  the 

45°  configuration  of  the  initial  locations  of  the  vortices.  (Note  that 
To  lies  in  the  second  quadrant  and  that  Ti  increases  positively  from 
To*)  Thus,  for  the  case  where  the  motion  is  periodic,  the  motion  in 
the  y,z  plane  is  completely  described  by  equations  (7),  (10),  (l4) 
and  (24)- 

If  the  motion  in  the  y,z  plane  is  not  periodic,  that  is,  if  G 
is  greater  than  4,  then  the  motion  of  vortex  1  is  described  by  equa¬ 
tion  (21).  Here  again  the  first  three  integrals  of  this  equation  are 
elliptic,  but  now  the  required  transformations  are 


and  it  is  recalled  that 


R  -  lyr  £ 

sh2  u.  =  sin2  cp1  =  - -i - 

1  1  rM2 


sn2  =  sin2  = 


r2-  U0- 1 


r2-q2 


Vs  -  Rg-Qa  _  M 

R2  G2 


(25) 
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It  should  he  noted  that  for  this  case  (aperiodic  motion)  when  T|2  is 
less  than  Q2,  all  four  integrals  become  imaginary.  That  is,  zx  is 


imaginary  if 


yi  _  i  |  <  q 
f  2  I  f 


This  condition  defines  the  asymptotes  of  the  vortex  paths  for  the 
aperiodic  motion.  This  can  he  seen  from  the  relative  paths  shown  in 
figure  1. 

If  the  integrals  of  equation  (21)  are  now  evaluated  for  this  case, 
the  resulting  expression  for  the  aperiodic  path  of  vortex  1  is 


—  =  -  -r— -  E(k,<Po)-  E(k,9x)  +  F(k,cp0)  -  F(k,9x)  + 

f  g2-i6  L  J  '  x 


2(G-4) 


^tan  <p0  ^-k2  sin2  cp0  -  tan  9i  */  1-k2  sin2  9i  ^  + 


1  £+b 

2  /  ~ 


tan  9o  -  tan  9i  J  +  ~ 


where  k  =  The  motion  is  now  completely  defined  by  equations  (7), 
(10),  (lb),  and  (27) . 

For  the  special  case  when  G  is  equal  to  k ,  all  the  integrals  are 
elementary  and  can  readily  he  evaluated  by  making  use  of  the  trigono¬ 
metric  substitutions 

f  f 

o 

cos  qPx  =  - cos  q>0  =  - -  (28) 

R  R 

The  final  solution  of  equation  (22)  is  then 


If.  \  1  / sin  9i  sin  9c 

r  sin  9,  -  sin  90  )+  -  (  - — 3 — - s-2- 

^  \  ‘v  2  \cosa  9i  cos2  90 


1  cos  9o  (l+  sin  9! ) _ 1_  f  sin  9i  sin  9o 

2  cos  9X  (l+  sin  9q  )  </2”  \  cos  9i  cos  9o 


2  V1  */J 
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Actually,  it  can  be  shown  that  this  result  is  obtainable  from  either 
equation  (24)  or  (27)  by  taking  the  limit  as  G  approaches  4.  It  is 
interesting  to  note  that  G  can  assume  only  positive  values  and  must 
lie  between  2  and  infinity  as  can  be  seen  from  its  definition  given 
with  equation  (l4) .  Thus,  the  vortex  paths  in  the  y, z  plane  have  been 
obtained  for  the  entire  range  of  the  parameter  G.  Figure  2  shows  the 
paths  of  vortices  1  and  2  for  several  values  of  G.  The  dashed  lines 
connect  corresponding  positions  of  the  two  vortices  at  successive  time 
intervals . 


Complete  motions  in  three  dimensions.-  In  order  to  complete  the 
three-dimensional  picture  of  the  four  vortex  filaments  trailing  behind 
the  cruciform  wings,  a  relationship  must  be  obtained  between  the  points 
of  the  paths  in  the  y,z  plane  and  the  distance  the  wing  has  traveled 
in  the  stream  direction.  This  can  be  done  by  the  relationship 


d 

« 


Ut 


=  u 


ay 


ay,  /at 


(30) 


Equations  (l4)  and  (l8)  can  again  be  used  to  express  v*  as  a  function 
of  y1  so  that  equation  (30)  can  be  integrated.  The  integral  thus 
becomes 


8jtGfU  , 

r 1  /  j 

2  2 

yi  (f-yi) 

r  ^ 
1  j 

r0  ^/^y1(f-y1)  )2 

[  ^^(f-yi)]5 

dyi 


(31) 


Again  using  the  substitution 


y,  =  n  +  a- 


equation  (31)  can  be  simplified  into  three  somewhat  different  forms, 
one  for  each  of  the  three  regimes  of  the  parameter  G. 


For  the  case  when  G  is  less  than  4  (periodic  motion) ,  equa' 
tion  (31)  can  be  written  as  2 


d  =  - 


8itfU 


7,-2 


rJ~G{G  +  4) 


J  f 

y0-2 


(f  -  .')■ 


dr] 


(Q2+T]2)  V  (R2-1]2)  (Q2+r|2) 


(32) 


*r 


» 


vhere 


R2  -  A 


4(g  +  4) 


and  d2  -  ^  -  G) 


4g 


tr 
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For  G  equal  to  4, 


With  the  aid  of  equations  (35),  (36),  and  (37),  the  positions  of 
the  four  rolled-up  vortices  are  completely  described3  for  all  distances 
behind  the  cruciform  wings.  Thus,  for  G  less  than  4,  equations  (7), 
(10),  (l4) ,  (24),  and  (35)  completely  determine  the  motion.  The  other 
cases  are  given  by  the  corresponding  equations. 

APPLICATION  TO  CRUCIFORM  TRIANGULAR  WING  AT  45°  BANK 

Determination  of  Initial  Vortex  Positions 
and  Vortex  Strengths 

The  equations  thus  far  developed  give  the  motions  of  the  rolled-up 
vortices  in  terms  of  their  initial  positions.  In  order  to  relate  these 
motions  to  physical  cases  involving  a  given  cruciform  wing  at  a 
specified  angle  of  attack,  the  initial  positions  of  the  four  vortices 
must  be  determined  as  a  function  of  the  lift  coefficient.  This  amounts 

3  "  "  ~  "  ' 

The  evaluation  of  the  circulation  T  will  be  discussed  in  the  next 

section. 
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to  the  determination  of  the  parameter  G  in  terms  of  the  physical 
quantities  such  as  lift  coefficient  and  aspect  ratio.  Therefore,  this 
section  will  be  concerned  with  the  initial  locations  of  the  four  vortices 
for  a  given  cruciform  wing  (banked  V50)  at  a  given  flight  condition. 

Since  the  vortices  trailing  behind  the  cruciform  wing  actually  leave 
the  trailing  edges  as  two  flat  vortex  sheets,  the  initial  locations  of 
the  four  rolled-up  vortices  are  somewhat  fictitious.  However,  if  each 
wing  of  the  cruciform  is  considered  as  a  plane  wing  in  sideslip,  the 
locations  of  the  rolled-up  vortices  from  each  wing  can  be  calculated 
from  the  span  ioading  curve.  Such  a  calculation,  then,  can  be  considered 
as  yielding  the  positions  of  the  four  vortices  of  the  cruciform  before 
any  interaction  has  taken  place  between  the  vortices  of  the  two  wings, 
namely,  immediately  behind  the  trailing  edges.  Hence,  in  order  to  calcu¬ 
late  the  initial  positions  of  the  four  vortices,  the  span  loadings  of  the 
separate  wings  must' be  determined. 

In  reference  9,  it  was  shown  that  the  load  distribution  on  each 
wing  of  a  pointed  low-aspect -ratio  cruciform  wing  with  no  body  is  given 
by  the  expression 


where  Ap  is  the  difference  in  pressure  between  corresponding  points 
on  the  upper  and  lower  surfaces  of  the  individual  wing. 

Equation  (38)  is  not  valid  for  angles  of  sideslip  greater  than  . 

dx 

The  span  loading  is  obtained  by  integrating  this  load  distribution  in 
the  chordwise  direction;  that  is 
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Equation  (39)  can  be  separated  into  two  integrals,  and  since  for  trian¬ 
gular  wings  ds/dx  is  a  constant,  the  resulting  span  loading  is 


clc  -  4a*  sQ 


^  +  d^fe i los 


1  + 


v^i-y2/ 


|y/so 


or,  since  ds/dx  =  ^  for  triangular  wings. 


c7c  =  4af  s 
l  o 


(40) 


(41) 


In  the  case  of  the  unyawed  plane  wing,  it  is  clear  from  symmetry 
considerations  that,  if  the  vortex  sheet  is  to  roll  up  into  two 
vortices,  all  the  vorticity  on  one  side  of  the  plane  of  symmetry  must 
eventually  become  rolled  up  into  one  vortex.  This  leads  to  the  conclu¬ 
sion  that  each  rolled-up  vortex  must  be  of  strength  equal  to  the  maximum 
circulation  which  in  this  case  is  the  circulation  in  the  plane  of  sym¬ 
metry.  A  corollary  of  this  conclusion  is  that  each  rolled-up  vortex 
contains  only  vortices  of  the  same  sense.  This  seems  to  be  in  accord¬ 
ance  with  existing  knowledge  of  the  behavior  of  vortex  sheets  as  well  as 
with  existing  theories  of  the  cumulation  of  vorticity  as  applied  to 
turbulence.  For  the  wing  in  sideslip,  whereas  there  is  no  symmetry,  the 
wing  can  again  be  considered  as  composed  of  two  segments,  each  having 
vorticity  in  only  one  direction  and  hence  each  producing  one  rolled-up 
vortex.  The  dividing  line  between  two  such  segments  is  at  the  spanwise 
station  of  maximum  circulation  so  that  such  a  division  leads  to  the 
conclusion  that  each  rolled-up  vortex  must  have  a  strength  equal  to  the 
maximum  circulation.  Thus,  the  vortex  strength  is  given  by 


^  ~  2  ^C^C^max  (^2) 

Further,  the  distance  between  the  vortices  is  given  by  the  fact  that  the 
lift  impulse  must  always  be  that  of  the  wing  itself.  Thus, 


2s r 


v  m  CL’S 

PU  r  "  (cZc)max 


(43) 


From  equation  ( 43 )  it  can  be  seen  that  the  distance  2s f  can  be 
represented  on  the  span  loading  curve  as  shown  in  the  sketch,  where  the 
area  under  the  rectangle  is  equal  to  the  area  under  the  curve. 
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As  seen  from  the  sketch,  the  actual  location  of  each  of  the  vortice 
still  requires  the  determination  of  the  distances  a*  and  a2  k  It  has 
already  teen  argued  that  all  the  vorticity  on  one  side  of  the  maximum 
circulation  must  eventually  he  contained  in  one  rolled-up  vortex.  Thus 
the  two  areas  can  be  equated  on  either  side  of  the  maximum;  that  is, 


(ai  +  ym)(czc): 


c7c  dy 


(aa  “  yin) (cZc)max  =  /  C2C  ^ 


(Note  that  ym  is  readily  calculated  by  successive  approximations.)” 

To  apply  the  above  information  to  the  cruciform  wing  at  45°  bank, 
one  merely  needs  to  note  that  for  this  case 

af  =  p  =  a  sin  45°  = 

v  2 


and  that  the  initial  position  of  vortex  number  1  in  the  banked  coordi¬ 


nate  system  is  now  given  by  y0  = 


It  is  further  pointed  out  that 


f  =  yi  +  72  = 


ax  +  a2 
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Since  the  lift  of  the  cruciform  wing  can  now  he  expressed  as 

L  =  pur  (ai  +  a  a)  =  |  pU2  CL  S  (1+5) 

and  the  wing  area  S  is  equal  to  4s02/A,  it  is  seen  that  the  circula¬ 
tion  T  is  given  by 

r  ..  JJ_ Sog  (46) 

U  ai  +  q.2  \  A  / 

Sufficient  information  is  now  at  hand  to  determine  the  initial 
positions  of  the  four  vortices  as  a  function  of  the  angle  of  attack  for 
a  particular  aspect  ratio  or,  more  generally,  as  a  function  of  a/A. 
However,  since  the  above  span  loadings  were  obtained  from  low-aspect- 
ratio  theory,  for  which  Cl  =  ~  Aa,  the  results  can  be  considered  as  a 

function  of  the  parameter  Cl/aA2.  It  can  be  shown  that  this  will  make 
the  results  more  general  in  that  they  will  now  apply  with  the  accuracy 
of  linearized  theory  rather  than  of  low-aspect -ratio  theory.  This  may  be 
seen  from  the  fact  that  the  loadings  in  linearized  theory  (see,  e.g., 
reference  12)  for  triangular  wings  are  simply  those  of  low-aspect-ratio 
theory  multiplied  by  a  constant  which  depends  upon  the  aspect  ratio  and 
the  Mach  number. 

The  theoretical  initial  positions  of  vortex  1  have  been  calculated 
and  are  plotted  against  C^/itA2  in  figure  3«  From  these  initial  posi¬ 
tions,  the  values  of  G  were  calculated  (equation  (l4)),  and  are 
plotted  in  figure  4.  The  limiting  angle  of  sideslip  of  equations  (38) 
and  (4l)  has  been  exceeded  somewhat,  as  indicated  by  the  dashed  portion 
of  the  curve  in  figure  3>  in  order  to  permit  some  interesting  observa¬ 
tions  regarding  the  indicated  trends.  The  constants  appearing  in  the 
equations  for  the  vortex  paths  have  now  been  completely  determined  and 
the  three-dimensional  paths  of  the  four  vortices  can  be  calculated  for 
any  lift  coefficient  and  aspect  ratio.  An  illustration  of  such  paths 
for  a  typical  case  is  presented  in  figure  5* 

If  one  recalls  the  types  of  motion  associated  with  the  various 
regimes  of  G,  figure  3  takes  on  added  significance.  The  passing  from 
periodic  motion  to  aperiodic  motion  simply  indicates  the  inward  move¬ 
ment  of  the  initial  position  of  the  upper  two  vortices  with  angle  of 
attack  until  their  velocity  in  the  z  direction  is  so  great  that  once 
they  pass  through  the  lower  vortices,  the  latter  never  catch  up  to 
complete  the  cycle.  If  this  inward  movement  were  to  continue  as  the 
angle  of  attack  increases  (as  shown  in  fig.  3)  then  at  =  0.244 

the  two  upper  vortices  would  have  coalesced,  leaving  only  the  two  lower 
vortices  which  would  then  travel  parallel  to  the  z  axis  as  in  the  case 
of  a  plane  wing. 
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Leapfrog  Distance 


As  an  indication  of  the  practical  importance  of  the  motion  of  the 
four  rolled-up  vortices  in  the  calculation  of  downwash,  a  distance 
behind  the  wing  which  is  characteristic  of  the  motion  will  be  calculated. 
This  distance  is  taken  as  the  distance  at  which  the  four  vortices  have 
become  collinear,  that  is,  the  distance  behind  Isjie  ;wing  at  which  the 
upper  two  vortices  are  just  passing  through  the  lower  two.  This  will  be 
referred  to  as  the  leapfrog  distance  dL  measured  from  the  wing  trailing 
edges. 

The  leapfrog  distance  can  be  calculated  directly  from  equation  (31) 
by  taking  the  upper  limit  of  integration  to  be  the  value  of  y1  in  the 
collinear  configuration.  This  value  of  yx  is  obtained  from  equa¬ 
tion  (l4)  by  setting  zx  2  equal  to  zero.  Solution  of  this  equation 
then  gives  the  limit  of  integration  yx  as  a  function  of  G  only  and 
therefore  of  CL/*A2.  The  definite  integral  of  equation  (31)  has  been 
evaluated  by  the  methods  discussed  previously  for  the  entire  range  of  G 
and  the  resulting  leapfrog  distances  are  plotted  (in  terms  of  chord 
lengths)  against  C^/flA2  in  figure  6.  For  purposes  of  comparison,  as 
well  as  to  provide  a  measure  of  the  usefulness  of  the  assumption  that 
the  four  vortices  are  fully  rolled  up  at  the  trailing  edge,  the  distance 
for  the  vortex  sheets  to  roll  up  (as  calculated  by  Kaden,  reference  1)  is 
also  shown  in  this  figure.  Since  the  distance  to  roll  up  as  shown  in 
figure  6  is  that  for  a  plane  wing,  it  should  be  taken  as  an  indication 
of  the  average  distance  for  the  four  vortex  sheets  to  roll  up.  Actually, 
because  of  the  asymmetric  span  loading,  the  lower  vortices  will  roll  up 
somewhat  faster,  while  those  from  the  upper  wing  panels  will  roll  up 
more  slowly. 

The  outstanding  point  to  be  noted  from  figure  6  is  that  the  phenom¬ 
enon  of  leapfrogging  can  occur  within  a  few  chord  lengths  at  reasonable  ■ 
lift  coefficients  for  low-aspect -ratio  triangular  cruciform  wings. 

The  calculations  are  expected  to  represent  the  physical  phenomenon  most 
accurately  when  the  leapfrog  distance  is  considerably  greater  than  the 
distance  for  the  vortex  sheets  to  roll  up.  The  curves  presented  in 
figure  6  have  been  terminated  at  a  value  of  Cp/itA2  of  0.176,  since 
above  that  value  (calculated  from  low-aspect -ratio  theory)  the  upper 
leading  edges  of  the  cruciform  wing  become  trailing  edges  and  the  span 
load  distributions  are  no  longer  given  by  equation  (^l) . 

Rolling  Up  of  the  Vortex  Sheets 


As  was  mentioned  at  the  outset  of  the  analysis,  the  assumption  that 
the  wings  of  the  cruciform  act  independently  of  one  another  is  not 
strictly  applicable  to  the  present  problem  because  once  the.  bound 
vortices  of  the  wings  become  free  (i.e.,  immediately  behind  the  trailing 
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edge)  the  two  vortex  sheets  influence  one  another  in  their  subsequent 
behavior.  Thus,  while  the  foregoing  analysis  has  treated  in  some  detail 
the  complete  motions  of  the  fully  rolled-up  vortices,  it  has  not 
considered  the  mutual  effects  of  the  vortex  sheets  during  the  rolling  — 
up  process  itself.  Although  a  complete  knowledge  of  the  rolling-up 
process  could  be  gained  only  by  a  detailed  numerical  analysis,  some 
important  questions  regarding  the  over-all  interaction  effects  can  be 
answered  by  examining  the  nature  of  the  vortex  sheets  in  their  initial 
cruciform  configuration.  Illustrated  schematically  in  the  sketch  is  the 

initial  vortex  config¬ 
uration,  showing  the 
span  loading  or  circu¬ 


lation  distribution, 
which  is  the  same  on 
each  wing,  and  the 
resulting  vorticity  with 
directions  of  rotation 
indicated  by  the  arrows. 
In  view  of  the  arguments 
presented  earlier,  it  is 
clear  from  the  sketch 
that  all  the  vortices 
contained  in  the  dis¬ 
tance  AO  will  roll  up 
into  a  single  vortex, 
and  that  all  those  in 
BC  will  roll  up  into 
another  single  vortex 
with  opposite  direction 
of  rotation.  Further¬ 
more,  the  strength  of 
the  latter  vortex  must 
clearly  be  equal  to  the 


maximum  circulation 


On  the  other  hand,  it  is  not  clear  just  what  will  become  of  the 
vorticity  contained  in  OB.  If  the  wings  did,  in  fact,  act  independ¬ 
ently,  then  OB  would  certainly  be  combined  with  AO  to  form  a  single 
vortex  of  strength  Lay.  However,  the  presence  of  a  plane  of  symmetry 
at  y  =  0  requires  that  AO  and  OB  now  be  considered  as  separate 
vortex  sheets  since  the  vorticity  at  0  vanishes  due  to  symmetry. 

In  fact,  closer  examination  of  the  portion  of  the  sheets  near  0  reveals 
that,  due  to  the  directions  of  rotation,  the  broken-line  sheet  AOB' 
must  move  away  from  the  sheet  A'OB,  thus  severing  AB  and  A'B’  at  0. 

It  is  doubtful,  then,  that  the  vorticity  originally  contained  in 
OB  and  OB'  will  ever  become  entrained  as  a  part  of  the  upper  two 
vortices.  It  will,  in  fact,  later  be  shown,  from  the  experiments  con¬ 
ducted  for  this  report,  that  this  amount  of  vorticity  actually  forms 
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a  separate  sheet  which  is  not  identified  as  part  of  any  of  the  four 
vortices  treated  thus  far. 


The  foregoing  observations  suggest  a  possible  improvement  on  the 
four-vortex  analysis  already  carried  out;  that  is,  the  strengths  and 
initial  positions  of  vortices  1  and  4  could  be  modified  to  take  into 
account  the  fact  that  the  vorticity  in  OB  and  OB'  actually  is  not 
contained  in  vortices  1  and  k.  The  amount  of  vorticity  in  OB 
(  ereafter  referred  to  as  rc)  and  the  resulting  modification  of  I\ 
is  readily  calculated  from  the  span  load  distribution,  as  shown  in  the 
sketch.  The  strengths  of  the  three  vortices  shown  are  simply  given  by 


ri  "  2  (cZc)o;  r2  =  \  (cZc)max 


rc  =  V1*  =  f  [(cZc)max  ~  (cZc)0  ~ 


,,  6  T°sltlons  of  "the  vortices  are  again  calculated  by  equating 

+ne^llft-lmPU^S\^ef°re  and  after  rolling  UV’  but  now  r-L  is  considered 
contain  only  the  vorticity  arising  from  its  side  of  the  wing;  that  is. 


ai!  (cZc)o  -  /  cic  dy 


ac  (CZC)0  +  (ym"ac)(c7c) 


c7c  dy 


(a2-ym)  (cZc)max  =  /  C2C  dy 
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Note  that  a2  is  unchanged  by  the  introduction  of  Tc  . 


With  the  aid  of  the  foregoing  sketch,  one  can  again  express  the 
cruciform-wing  lift  in  terms  of  the  initial  vortex  positions 


L  =  V2  pU  ri  (aj.*  +  ac)  +  r2  ( 


l.i 
j  2 


pU  CLS  (49) 


and  conclude  that 


JT  Sn2  fa 


u  [y(a>i'  +  ac)  +  a2  —  ac  ] 


where 


The  original  four -vortex  prob- 
lem  has  actually  now  been  replaced 
^  by  the  six -vortex  problem  illus¬ 

trated  in  the  sketch.  However,  due 
^4  1  to  the  relatively  small  strength  of 

as  veH  as  to  the  close  prox¬ 
imity  of  the  two  opposing  vortices 
*  Fc,  the  influence  of  these  addi- 

- y  tional  vortices  on  the  motions  of 

1  c  7>i  rC  the  other  four  is  expected  to  be 

V  small.  An  iterative  approach  to 

the  exact  solution  can  be  obtained 
by  neglecting  this  influence  once 
„  the  initial  position  and  strength 
"3  V  -J  of  r±  have  been  modified  to  allow 

for  Tc.  Thus  the  problem  is  again 
simplified  to  that  of  the  motion  of 
four  vortices,  but  now  they  have 
unequal  strengths.  (The  ratio  of 
the  strengths  Vx  to  F2  is  plotted  against  Cp/itA2  in  fig.  7-) 

This  fact  alone  complicates  the  mathematics  to  the  point  that  the 
integrations  can  no  longer  be  carried  out  except  by  numerical  methods. 
Therefore,  the  complete  motions  will  not  be  determined,  but  the  leapfrog 
distance  will  be  recalculated  in  order  to  illustrate  the  effect  of  the 
new  assumptions  on  the  results.  The  leapfrog  distance  is  readily 
checked  experimentally*  as  will  be  shown  in  the  experimental  section. 
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Modified  Four-Vortex  Calculation 


The  equations  for  the  relative  motion  of  the  four  vortices  can 
again  be  set  up  by  the  use  of  the  general  equations  (3)  and  (4)  with  a 
new  statement  of  the  center-of -gravity  rule  in  which  the  vortex 
strengths  are  no  longer  taken  to  be  equal.  That  is, 


yc.g. 


£yiri  riyi  +  rayg 

sri  ri  +  r2 


e  =  constant 


(51) 


If  it  is  noted  that  in  the  initial  configuration 

zi} 2=  yi  +  Yi  =  y0;  Ys  =  e  (i  +  7)  -  ryQ  (52) 


then  equations  (3)  and  (4)  can  be  used  as  in  the  original  analysis, 
resulting  in  the  following  expression  (in  dimensionless  form)  for  the 
relative  vortex  paths: 


The  velocity  v-l  is  recalculated  in  the  same  manner  as  before,  with 
the  result 
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( 5<0 

If  this  egression  is  now  used  in  equation  (30)  to  determine  the 
downstream  distance  d,  the  resulting  (dimensionless)  form  of  the 
integral  is 


It  is  readily  ?seen,  without  substituting  the  egression  for  zx  2/e 
as  given  by  equation  (53) >  that  analytical  evaluation  of  the  integral 
of  equation  (55)  would  be  extremely  difficult.  This  integration  was 
carried  out  numerically  with  the  definite  upper  limit  as  given  by  the 
collinear  configuration  of  the  four  vortices  (i.e.,  zx  2  =  0)  and  the 
resulting  values  of  the  leapfrog  distance  d^  are  plotted  on 
figure  6.  It  is  noted  that  the  resulting  curve  is  higher  than  that 
calculated  in  the  original  analysis.  This  is  due  primarily  to  the 
smaller  strengths  of  the  upper  two  vortices  and  to  the  consequently 
longer  time  required  for  them  to  pass  through  the  lower  two  because 
of  reduced  induced  velocities  of  one  upon  the  other. 


WATER-TANK  EXPERIMENTS 

\ . 

* 

The  motions  of  the  vortex  sheets  behind  a  triangular  cruciform 
wing  of  aspect  ratio  2  were  observed  experimentally  by  means  of  a 
water  tank  (fig.  8) .  The  model  was  mounted  on  a  vertical  track  which 
was  driven  at  uniform  speed  into  the  tank,  while  photographs  of  the 


■water  surface  were  recorded  “by  a  motion  picture  camera.  The  traces  of 
the  vortex  sheets  were  made  visible  by  applying  fine  aluminum  powder  to 
the  trailing  edges  of  the  model.  The  model  had  a  span  of  8  inches  and 
was  made  of  0.050-inch  sheet  metal. 

In  general,  the  purpose  of  the  water-tank  studies  was  simply  to 
illustrate  the  motions  analyzed  in  the  earlier  portions  of  this  report. 
(A  typical  series  of  photographs  enlarged  from  the  moving  pictures  is 
presented  In  fig.  9*)  However,  it  was  a  simple  matter  to  obtain 
experimental  values  of  the  leapfrog  distance  by  means  of  a  tape  which 
moved  with  the  model  track  and  recorded  on  the  movie  film  the  distance 
traveled  by  the  wing.  The  results  of  such  experiments  are  presented  in 
figure  6  for  several  angles  of  attack  and  it  is  observed  that  the  agree¬ 
ment  with  the  modified  four -vortex  calculation  is  remarkably  good. 

The  formation  of  a  separate  vortex  sheet  arising  from  the  center 
portion  of  the  cruciform  wing,  as  discussed  in  a  previous  section,  can 
be  seen  in  the  photographs  of  figure  9-  This  is  the  phenomenon  that 
has  been  taken  into  account  in  the  modified  theory. 

It  is  interesting  to  note  (fig. 9)  that  the  two  lower  vortices 
extend  downstream  in  nearly  the  free-stream  direction  as  indicated  by 
their  positions  relative  to  the  wing -tip  markers.  This  was  also  true 
for  the  calculated  paths  of  figure  5* 

An  important  feature  of  the  actual  flow  field  as,  distinguished  from 
the  simplified  model  used  for  the  analysis  is  the  persistence  of  the 
vortex  sheets  between  the  rolled-up  vortices  as  seen  in  figure  9* 

Because  of  this  fact,  the  theoretical  paths  of  the  vortices  are  not 
expected  to  be  accurate  at  distances  behind  the  wing  greater  than  dL' 
since  vortices  1  and  k  will  then  begin  to  become  entrained  in  the  outer 
windings  of.  the  sheet  that  constitutes  vortices  2  and  3*  However,  the 
magnitude  of  this  effect  will  depend  upon  the'  rate  of  rolling *up  of 
the  vortex  sheets  and  it  would  be  expected  to  be  most  serious  at  the 
lower  lift  coefficients  (or  higher  aspect  ratios)  where  the  motions 
being  considered  are  not  of  practical  concern. 


CONCLUDING  REMARKS 


A  detailed  analysis  has  been  made  of  the  motions  of  the  four  fully 
rolled-up  vortices  trailing  behind  a  cruciform  wing  which  is  banked  45° • 
Equations  have  been  developed  describing  the  paths  of  the  four  vortices 
in  three  dimensions,  and  calculations  have  been  made  of  the  distance 
behind  the  wing  at  which  the  upper  two  vortices  leapfrog  through  the 
lower  two.  The  latter  calculations  were  confined  to  wings  of  triangular 
plan  form.  The  simplified  analysis  presented  in  the  early  portions  of 
the  report  has  been  modified  to  account  for  the  fact  that  all  the 
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vorticity  shed  from  the  cruciform  wing  is  not  actually  contained  in  the 
four  rolled-up  vortices. 

It  was  found  that  the  leapfrog  distance  decreased  with  increasing 
lift  coefficient  and  increased  with  increasing  aspect  ratio.  For  angles 
of  attack  up  to  a  certain  critical  value  of  CL/jtA2,  the  vortex  motion  is 
periodic  with  downstream  distance,  while  above  the  critical  value  the 
motion  is  aperiodic. 

From  the  fact  that  the  leapfrog  phenomenon  can  occur  within  two 
chord  lengths  of  the  wing  trailing  edges,  it  is  clear  that  downwash 
calculations  even  at  distances  less  than  one  chord  length  behind 
low-aspect-ratio  cruciform  wings  must  take  into  account  the  vortex 
motions  considered  in  this  report.  Once  the  positions  of  the  vortices 
are  known,  there  are  methods  available  for  approximating  the  downwash, 
and  corrections  can  be  made  for  the  viscous  vortex  cores  by  assuming 
that  they  rotate  as  solid  bodies. 

For  angles  of  bank  other  than  45°,  it  is  doubtful,  due  to  the  lack 
of  symmetry,  that  calculations  of  the  type  presented  here  could  be  made 
by  other  than  numerical  procedures. 


Ames  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Moffett  Field,  Calif.,  Oct.  5,  1951. 
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Vortex  paths  in  the  yz  plane. 
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,1  paths  of  four  vortices  .trailing  behind  cruciform  wing  banked 
aspect  ratio  2.  Cl  =  1.^1. 
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Figure  7.-  Variation  of  relative  vortex  strengths 

with  lift  coefficient. 
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Figure  8.-  Water  tank  with  cruciform-wing  model. 
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d/c  =  0.33  d/c  =  1.38 


Figure  9*-  Photographs  of  wake  at  various  stations  behind  a  triangular 
cruciform  wing  of  aspect  ratio  2.  Cl  r;  0.66. 
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a/c  =  2.83 


a/c  =  4.26 


(b)  Stations  7  to  11. 


Figure  9--  Concluded. 
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